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Abstract 

A set of 55 clarinet reeds is observed by holography, collecting 2 series of measurements 
made under 2 different moisture contents, from which the resonance frequencies of the 15 
first modes are deduced. A statistical analysis of the results reveals good correlations, but 
also significant differences between both series. Within a given series, flexural modes are not 
strongly correlated. A Principal Component Analysis (PCA) shows that the measurements 
of each series can be described with 3 factors capturing more than 90% of the variance: 
the first is linked with transverse modes, the second with flexural modes of high order and 
the third with the first flexural mode. A forth factor is necessary to take into account 
the individual sensitivity to moisture content. Numerical 3D simulations are conducted by 
Finite Element Method, based on a given reed shape and an orthotropic model. A sensitivity 
analysis revels that, besides the density, the theoretical frequencies depend mainly on 2 
parameters: El and Glt- An approximate analytical formula is proposed to calculate 
the resonance frequencies as a function of these 2 parameters. The discrepancy between 
the observed frequencies and those calculated with the analytical formula suggests that 
the elastic moduli of the measured reeds are frequency dependent. A viscoelastic model is 
then developed, whose parameters are computed as a linear combination from 4 orthogonal 
components, using a standard least squares fitting procedure and leading to an objective 
characterization of the material properties of the cane Arundo donax. 

1 Introduction 

Clarinettists experience every day the crucial importance of clarinet reeds for the quality of 
sound. Their characterization is a real challenge for musicians who wish to obtain reeds that are 
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suited to their personal needs. The present paper address this complex field of research, but its 
scope is restricted to the development of an objective method for a mechanical characterization of 
single reeds of clarinet type. >From the shape and the resonance frequencies of each individual 
reed (measured with heterodyne holography), we intend to deduce the mechanical properties 
of the material composing it. A subsequent study should then examine how these mechanical 
properties are correlated with the musical properties of the reeds. 

Natural materials, as wood or cane, are often orthotropic and exhibit a different stiffness 
along the grain (longitudinally) as in the others directions. The problem is then obviously mul- 
tidimensional. Nevertheless, reed makers classify their reeds by a single parameter: the nominal 
reed "strength" (also called "hardness"), in general from 1 to 5, which basically reflects the stiff- 
ness of the material (cane, Arundo donax L.), since all reeds of the same model have theoretically 
the same shape. The method of measurement is generally not publicized by manufacturers, but 
this "strength" is probably related to the static Young modulus in the longitudinal direction El- 

"Static" (i.e. low frequency) measurements of the elastic parameters of cane are available 
in the literature, for instance Spatz et al. [Tj. A viscoelastic behavior has been reported in 
experimental situations (see e.g. Marandas et al. [2], Ollivier [3] or Dalmont et al. [I]) and this 
fact seems generally well accepted in wood sciences and biomechanics (for instance Speck et al. 
[SJ [5] ) - Marandas et al. proposed a viscoplastic model of the wet reed. Viscoelastic behavior for 
cane was already demonstrated by Chevaux [7], Obataya et al. [5J H2 [TU1 [H] and Lord [12]. These 
authors study only the viscoelasticity of the longitudinal Young modulus El, leaving aside the 
fact that of the shear modulus in the longitudinal/tangential plane Glt- Furthermore, they give 
no really representative statistics about the variability of the measured parameters. 

Different authors (among them Casadonte [15 1 114 ] , Facchinetti et al. [T51IT5] and Guimezanes 
|17j ) modeled the clarinet reed by Finite Elements Method (FEM) and computed the first few 
eigenmodes. They chose appropriated values of the elastic parameters in the literature, ignor- 
ing however viscoelastic behavior. The goodness of fit between observations and model was of 
secondary importance, except for Guimezanes. 

This latter author built a 2-D elastic model of the reed with longitudinally varying parameters. 
He fitted his model quite adequately with his observations (only 5 resonances were measured), 
but the fitted parameters seem not really plausible physically. His model didn't respect the 
assumption of a radial monotonically decrease of stiffness from the outer side to the inner side 
of the cane. Under such conditions, the frequency of the first resonance would increase in 
comparison to homogeneous material, and not decreased, as observed experimentally. 

The observation of mechanical resonance frequencies can be achieved by different methods. 
The methods used by Chevaux, Obataya and Lord are destructive for the reed, which cannot 
be used for further musical tests. On the contrary, holography is a convenient non-destructive 
method, the reed being excited by a loudspeaker. For instance Pinard et al. [THj measured with 
this method the frequency of the 4 lowest resonances and focused their attention on the musical 
properties of the reeds. 

The digital Fresnel holography method was used by Picart et al. 19, 20 and Mounier et al. 
[2"T] to measure high amplitude motion of a reed blown by an artificial mouth. Guimezanes [T7] 
used a scanning vibrometer. 

Recent technological developments provide very efficient and convenient measurements with 
holography, without having to manually identify the modes of resonance and to be satisfied with 
a single picture of their vibration: in a few minutes hundreds of holograms are acquired showing 
the response of a reed for many frequencies. The temperature and the moisture content can be 
considered as constant during a measurement serie^l The Sideband Digital Holography technique 

1 The significantly lower correlations between resonance frequencies (compared to our data) shows that it was 
probably not the case in Pinard's study. This fact may also reflect an unprecise determination of the resonance 
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provides additional facilities (see I2.1.ip . 

Generally, the physicist chooses a model in order to validate it by observations. In the present 
study, the complexity of the problematic forced us to adopt the reverse attitude: We observe 
the mechanical behavior of clarinet reeds with a statistically representative sample and exploit 
afterward the statistical results for establishing a satisfactory mechanical model designed with a 
minimal number of parameters. 

In section [5] the measurement method is presented. The experimental setup is described in 
section 12.11 and the method for observing resonance frequencies is detailed in 12.21 The results 
for 55 reeds are given in section [3] (statistics, correlations and Principal Component Analysis 
(PCA)[22). 

In section the development and the selection of a satisfactory mechanical model with 
minimal structure is described. First, a numerical analysis of the resonance frequencies of a reed 
assumed to be perfectly elastic is done by Finite Element Method (FEM), using the software 
Catia by Dassault systems, and a metamodel computing the resonance frequencies from elastic 
parameters is given in section l4TTl This allows to solve the inverse problem in a fast way. However, 
because the elastic model is not very satisfactory, viscoelasticity has to be introduced and some 
parameters are added to the model in section FOl The viscoelastic model has however too many 
degrees of freedom, according to PCA. Consequently, the viscoelastic parameters of the model are 
assumed to be correlated and PCA indicates that these parameters can be probably reconstructed 
from 4 orthogonal components, as a linear combination, by multiple regression fsection l4.3[) . The 
relationships between the components and the viscoelastic parameters is given, and finally the 
resulting values for these parameters are discussed in section 14.41 and compared with the results 
of the literature. 

2 Observations by Sideband Digital Holography 

2.1 Experimental setup 
2.1.1 Holographic setup 

The experimental setup is shown schematically in Fig. [T] A laser beam, with wavelength 
A = 650 nm (angular frequency uil) is split into a local oscillator beam (optical field E^o) and 
an illumination beam (Ei); their angular frequencies u>lo and wj are tuned by using two acousto- 
optic modulators (Bragg cells with a selection of the first order diffraction beam) AOM1 and 
AOM2: lulo = + ^AOMX an d Wj = lul + ^AOM2, where u>aomi,2 — 2ir x 80 MHz. The first 
beam (LO) is directed via a beam expander onto a CCD camera, while the second beam (I) is 
expanded over the surface of the reed, which vibrates at frequency / . The light reflected by the 
reed (field E) is directed toward the CCD camera in order to interfere with the LO beam (Elo)- 
4 phases were used (phase shifting digital holography) and we select the first sideband of the 
vibrating reed reflected light by adjusting uaom 1,2 to fulfil the condition: ujaomi — ^>AOM2 — 
+ /ccd)/4, where fcCD is the CCD camera frame frequency. The complex hologram 
signal H provided by each pixel of the camera, which is proportional to the sideband frequency 
component of local complex field E, is obtained by 4-phases demodulation: H = (Iq — I2) +j(h — 
I3) where Iq . . . I3 are 4 consecutive intensity images digitally recorded by the CCD camera, and 
j 2 — — 1. From the complex hologram H, images of the reed vibration are reconstructed by a 
standard Fourier holographic reconstruction calculation [23]. These holographic reconstructed 
images exhibit bright and dark interference fringes. Counting these fringes provides the amplitude 

frequencies. 
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Figure 1: Holographic setup. L: main laser; A0M1, AOM2: acousto-optic modulators; M: 
mirror; BS: beam splitter; BE: beam expander; CCD: CCD camera; LS: loudspeaker exciting 
the clarinet reed through the bore of a clarinet mouthpiece at frequency / = lo/2tt. 

of vibration of the object (in the direction of the beam) , which depends on the wavelength A of 
the laser, and on the first Bessel function J\, for instance ±95nm for the first, ±770nm for the 
5 th and ±1.6^m for the 10 th maximum (bright fringes) [2~4"l I25pl . 
This method has 3 main advantages: 

(i) The time for data acquisition is very short, about 3 minutes for recording 184 holograms, 
including holographic reconstruction. 

(ii) The signal to noise ratio is significantly better than with traditional technology, particu- 
larly through the elimination of signal at zero frequency. 

(iii) The visualization of large-amplitude vibration (order of magnitude: 0.1 mm) is possible 
by using high harmonics orders (up to several hundred times the excitation frequency). 

2.1.2 Reed excitation 

The reed was excited by a tweeter loudspeaker screwed onto an aluminium plate, connected to a 
clarinet mouthpiece. The lay of this mouthpiece was modified to be strictly flat. A plastic wedge 
of uniform thickness has been inserted between the lay and the reed, longitudinally to the same 
height as the ligature ( Vandoren Optimum) , allowing free vibrations of the entire vamp (length: 
about 38 mm), see Fig. [T] This ensures precise boundary conditions, avoiding any dependence 
to deformations of the reed. The repeatability of the longitudinal placing of the wedge and of 
the reed was ensured by a Claripatch ring [21)] . 
This setup requires some comments: 

(i) The reed is excited exclusively through the bore of the mouthpiece. 

(ii) The pressure field in the chamber of the mouthpiece was not measured. Like for a real 
instrument, the edges of the reed (protected by the walls of the chamber) are subject to a pressure 
field, which is probably lower than the pressure acting on the rest of the vamp. 

2 The original notation from the cited paper is kept. This notation is only valid for this paragraph. 



4 



Figure 2: Typical holographic patterns of series A (asymmetrical sinusoidal excitation: the left 
side of the reed is more strongly excited than the right side). Frequency range: 1.4 to 20 kHz by 
steps of 25 cents (181 pictures ordered from left to right, continued on the next row; the tip of the 
reed is down on each picture). Some modes are easily identified: FX (1st row, 19th picture), Tl 
(3rd row, 10th picture), F2 (4th row, 7th picture), T2 (4th row, penultimate picture), XI (5th 
row, 5th picture) etc ... Modes T3 and F3 are almost at the same frequency (6th row, 4th and 
5th pictures, probably). The last picture of the 3rd row corresponds to an acoustic resonance 
of the excitation device. It is present on all holograms of both series at the same frequency 
(examine Fig[3]). The excitation amplitude exponentially increases until the 73rd picture (3rd 
line, 21th picture), being hold constant beyond. 

(iii) The boundary conditions are very different from those of a real instrument (no curved 
lay, no contact with the lip). In addition, the reed was not moistened for the measurement. 

(iv) The excitation device is almost closed. The acoustical resonances of the excitation device 
are unknown, but may quite easily be deduced by comparing different measurements, because 
they are always present at the same frequency. 

2.2 Observation of resonance frequencies 
2.2.1 Experiment protocol 

55 clarinet reeds of model Vandoren V12 were purchased in a music shop: 12, 12, 20 and 11 
reeds of nominal forces 3, 3^, 4 and 4^, respectively. 29 reeds were used for two preliminary 
studies in order to develop the measurement protocol. Each of these reeds was played a total 
of some tens of minutes, spread over several weeks before measurement with the final protocol. 
The other 26 reeds were strictly new by measurement, which was performed immediately after 
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Figure 3: Typical holographic patterns of scries B (symmetrical sinusoidal excitation). The 
modes XI and T3 cannot be distinguished anymore. Notice that IT is less marked than un- 
der asymmetrical excitation (for some reeds even difficult to identify) and that the pattern has 
a significant flexural component, strongly dependent of the lateral placing of the reed on the 
mouthpiece. Notice that the symmetry of the patterns near Tl depends on the excitation fre- 
quency 



G 




Figure 4: (a) Flexural modes FX, F2, F3, FA and Fh. Side view, (b) transversal modes Tl, 
T2, T3, T4, T5. Front view, (c) Generic modes: XI, X2, X3, X4, X5 and X6. View from 
above. The intersections of nodal lines with the sides of the reed are symbolized by the blue 
dots. 

package opening (for 21 of them with the new hermetically scaled package by Vandoren, ensuring 
a relative humidity between 45 and 70%, according to the manufacturer), without moistening 
the reed. 

Each reed was subject to 2 series of measurements: 

• Series A (asymmetrical excitation: see FigJ2J): the right half of the mouthpiece chamber 
was filled with modelling clay to ensure a good excitation of antisymmetrical modes. 184 
holograms were made ranging from 1.4 to 20 kHz (sinusoidal signal), by steps of 25 cents. 
The amplitude of the excitation signal was exponentially increased in the range 1.4 to 4 
kHz, from 0.5 to 16 V, then kept constant at 16 V up to 20 kHz. This crescendo limits the 
amplitude of vibration of the first two resonances of the reed. The temperature was not 
measured (about 20°C). 

• Series B (symmetrical excitation: see FigJ3J): the modelling clay was removed. The 
protocol is otherwise identical to this of the first series. The reeds were inadvertently 
exposed during one night to the very dry and warm air from the optical laboratory between 
the two series of measurements. The reeds lost between 2 and 4% of their mass. In what 
follows we try to interpret the influence of this fact. The temperature was around 23-25°C. 

2.2.2 Nomenclature of normal modes 

Distinguishing 3 morphological classes, we classify the modes of a clarinet reed as follow : i) The 
"flexural" (or "bending", or "longitudinal") modes, listed below F, whose frequencies mainly 
depend on the longitudinal Young modulus (El) and polarized mainly in the z axis, ii) the 
"transversal" (or " torsional", or "twisting" ) modes, listed below T, mainly dependent on the 
shear modulus in the longitudinal / tangential plane (Glt), and iii) the "generic" (or "mixed") 
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Figure 5: The clarinet reed: coordinates system, (a) top view, (b) front view, (c) side view. 
x, y, z: Cartesian axes of the object. L, T, R: axes of the orthotropic material (L: longitudinal, 
T: tangential, R: radial). In the software Catia used for the simulations, the orthotropic model 
is Cartesian and not cylindrical. Therefore an exact equivalence between x, y, z, and L, T, R 
respectively can be assumed since we observed no important deviation between the direction of 
the grain and the axe of symmetry of the reed. 



modes, listed below X, sensitive to both moduli El and Glt (see Fig0]). A subclass of fiexural 
modes may be distinguished: the "lateral" modes (listed below L), polarized mainly in the y axis 
(see Fig. [5]). These modes were not observed in our study. 

The modes have been numbered after the order of increasing frequencies from a preliminary 
modal analysis we performed. In our analysis, however, the identification of a mode is based 
upon morphological criteria. As a matter of fact, the mode number and the order of observed 
frequencies are not necessarily identical for all reeds. 

Strictly speaking, the optical method only allows to observe the resonance frequencies of the 
reed and not the eigenfrequencies. Therefore the observed deformation patterns are a priori not 
identical to the eigenmodes of the reed. Nevertheless in practice no major differences have be 
found between the computed eigenmodes (see section 2]) and the observed or computed defor- 
mation for a forced asymmetrical excitation at the corresponding frequency. For this reason we 
use the terminology "mode" for the maximum amplitude of the response of the reed to a forced 
excitation. This is somewhat abusive, because the small shift between the resonance frequencies 
due to damping and the eigenfrequencies computed by FEM, without damping, is ignored. Be- 
sides damping, the acoustic load is also able to shift the resonance frequencies. We assume that 
this discrepancy is approximatively the same for all reed. 

2.2.3 Analysis of holograms; mode identification 

More than 30000 holograms were made for this study and analyzed as follows: The picture 
where the number of interference fringes is locally maximum is determined. For some cases, 
we chose the hologram that is most similar to our numerical simulations (by Finite Element 
Method, see section 0]) or to other holograms (see FigJS]) . The holograms corresponding to an 
acoustical resonance of the system, present at the same frequency (4309 Hz) for all reeds, have 
been eliminated. The identification of the different patterns to those calculated by FEM was 
often quite simple. An exception have been encountered for F"& and T3, whose frequencies were 
often so close that our identification is sometimes uncertain. More sophisticated techniques 
would certainly solve this problem. Notice that other boundary conditions (e.g. with clamping 
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Figure 6: Qualitative comparison with FEM computation, and typical variability of the experi- 
mental results, (a) First row: Quasi-static pattern (at 605 Hz, with strong excitation by the LS). 
2nd to 5th row: Flexural modes FX, F2, F3 and FA. Leftmost column: numerical simulation of 
eigenmodes by FEM. Columns 2 to 7: Arbitrary selection representing the observed variability 
(The first two rows correspond to the same selection of reeds). Notice the marked asymmetries 
and the differences in the curvature of the interference fringes near the tip of the reed, (b) Rows 
1 to 5: transversal modes TX, T2, T3, T4 and T5 (probably). Columns: see (a), (c) Rows 1 
to 6: Generic modes XX, X2, X3, XA, X5 and X6. The identification with X6 is sometimes 
unlikely. Columns: see (a). 



closer to the tip of the reed) would also easily separate these two modes. The frequency of some 
higher modes could not always be measured, either because their frequency was beyond 20 kHz, 
or because their pattern could not be clearly identified. 



3 Statistical analysis of resonance frequencies 

4 flexural, 5 transversal, and 6 generic modes have been identified, namely all 15 first modes of 
the reed, excluding lateral modes. This number is significant, compared to the 4 modes detected 
by Pinard et al. [TH]- The 6th mode (LI) could not be identified, as it is a lateral mode (flexural 
mode moving mainly in the y axis), not excited by our loudspeaker. We tried to observe it by 
rotating the mouthpiece to the side, without success. Notice that higher modes could probably 
be identified using an ultrasonic loudspeaker. 

3.1 Statistics 

The statistics are displayed on Table [T] For 14 measurements of resonance frequencies of the two 
series, identification of the total number of reeds (55) have been done. For other measurements, 
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Table 1: Observed resonance frequencies, sorted by frequency (in Hz). Na (resp. iVe): Number 
of identified pattern for each mode of series A (resp. B). /i: Mean value of the resonance 
frequency a: Standard deviation, min: Minimum, max: Maximum. Total: total number of 
identified patterns for each series. 



identification has been done only for a part of this number. The value of the ratio of the standard 
deviation a to the mean value fx, i.e. the relative standard deviation, is found to be between 
2 and 5% (about 1/3 tone). If we admit Gaussian distribution for the measured frequencies, 
99% of the observations typically range about ±1 tone (±200 cents) around the mean value (i.e. 
[i ± 3<t), for all frequencies. 

The identification of the mode X6 is uncertain: it seems to appear for frequencies lower than 
those of our simulations. Mode T5 is on the limit of the range we studied: this explains the 
small value of the standard deviation. 

Between series A and B, the flexural modes Fl to FA lower their midrange, while the transver- 
sal modes slightly increase it. The difference between the two series probably lies mainly in the 
drying of the reeds, and this seems to have a statistically significant effect. This is surprising, 
because drying decreases the density of the reed, and theoretically this should proportionally 
increase all frequencies. In addition, according to Obataya et al. [11] . drying is expected to 
increase E' L (at least around 400 Hz), which should also increase the resonance frequencies. 
However Chevaux [7] observed that drying diminishes E' L for material extracted from the inner 
side of the cane and augments slightly E' L for material extracted nearer from the outer side (for 
cane suitable for oboe reeds), at least in the frequency range 100-500 Hz. 

The hypothesis of an influence of the excitation method on the resonance frequencies seems 
unlikely, as well as the hypotheses of a poor reproducibility of the position of the reed on the 
mouthpiece between measurements or of the modification of the acoustic load, due to the mod- 
elling clay. 
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3.2 Correlations 



Linear correlations have been computed between all possible couples of variables in a usual way 
(all the variables from Table [TJ except X5 and T5 of series A, i.e. 13 variables for series A and 
15 for series B, plus the nominal reed strength). Results of series A for the mode Fl are denoted 
AF1, and similarly for the other results. The following 12 pairs have a correlation greater than 
0.9: AFl/BFlU, AT1/BT1, AT2/BT2, AT3/AT4, AT3/BT3, AT3/BT4, BT3/BT4, AT4/BT3, 
BX1/BX3, AT4/BT4, BT4/BX6 and BX3/BX4. 54 pairs of variables have a correlation between 
0.8 and 0.9, 50 other pairs between 0.7 and 0.8 and 262 other pairs, below 0.7. 

Between the two series, the correlation is excellent for corresponding transversal modes (Tl 
(i.e. AT1/BT1): 0.97, T2: 0.97, T3: 0.96 and T4: 0.98), and generally good for corresponding 
generic modes (XI: 0.87, X2: 0.84, A3: 0.87 and A4: 0.55). For flexural modes, the correlation 
is good for Fl and progressively lower for increasing mode order (Fl: 0.92, F2: 0.66, F3: 0.57 
and F4: 0.47). 

Within the same series, on the contrary, there is a poor correlation between AF1 and all 
measurements of series A, and similarly for BF1 and series B. This is striking: the two best 
correlated variables are AF2 and AX1 (0.73 and 0.49, respectively) for AF1 and BF2 and BT1 
(0.63 and 0.49, respectively) for BF1. Moreover, these correlations are quite low among all 
flexural modes: see Table [5] This fact is discussed in t)4.4l 

The nominal reed strength correlates at 0.7 with ATI and AX4. We expected a better 
correlation with Fl (only 0.6). This is surprising, since the reeds were probably sorted by a 
quasi-static bending method by the manufacturer. This would mean that the storage modulus of 
El at very low frequency is not well correlated with its value at the frequencies of the measured 
resonances. The influence of density has also to be considered. However, Obataya et al. [§] 
observed a good correlation between density and El- This point has to be investigated (see also 

3.3 Principal component analysis 

Principal Component Analysis (PCA) is mathematically defined as an orthogonal linear transfor- 
mation transforming the data to a new coordinate system, such that the greatest variance by any 
projection of the data comes to lie on the first coordinate (called the first principal component or 
first factor), the second greatest variance on the second coordinate, and so on [27 . Theoretically 
PCA is the optimum linear transform for given data in terms of least squares. PCA is based 
upon the calculation of the eigenvalue decomposition of the covariance (or of the correlation) 
matrix (see e.g. [22]). 

A PCA has been performed using the FACTOR module of SYSTAT [28 . The 14 variables 
(observed frequencies) presenting complete measurements for all reeds have been selected (all 
variables having 55 identified pattern, Na or Nb = 55, see Tabled]). Frequencies are rated in 
cents. 

The 4 largest eigenvalues have been selected. They capture 91.2% of the total variance of 
our sample (respectively 53.6%, 21.4%, 10.8% and 5.4% for each factor). A fifth factor would 
capture only 2.5% variance more. The 14-dimensional data have been linearly projected onto a 
4-dimensional factor space. 

The factor space can afterwards be orthogonally rotated, for instance for maximizing the 
correlations between rotated factors and observed variables. In the studied case, no a priori 
knowledge about the orientation of the factor space is available. For an easy comparison, using 

3 AF1/BF1 means AF1 versus BF1. The correlations are computed between AFl[n] and BFl[n], for n = 1 to 
N, where N = 55; missing observations are deleted. 
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0.177 


AX3 


0.761 


0.472 


0.090 


0.217 


AX1 


0.740 


0.451 


0.356 


0.115 


AF3 


0.223 


0.891 


076 


080 


AF2 


0.143 


0.791 


0.519 


0.064 


AF1 


0.085 


0.360 


0.870 


-0.141 


BF1 


0.104 


0.385 


0.835 


0.209 


BF3 


0.059 


0.596 


0.050 


0.755 


BF2 


0.147 


0.561 


0.290 


0.710 


ei 


0.958 


-0.080 


0.105 


0.005 


e2 


0.059 


0.969 


0.082 


0.098 


e3 


-0.075 


-0.073 


0.971 


0.042 


e 4 


-0.009 


-0.081 


-0.039 


0.979 



Table 2: Correlation (loadings) between rotated factors from PCA and : i) variables (mea- 
sured resonance frequencies) or ii) components from viscoelastic model (e[n], see § I4.3.1|) . for 
comparison, sorted in reverse order of magnitude. In bold: greater correlation for each variable. 

the VARIMAX algorithm, we choose to maximize the correlations between rotated factors and all 
available variables (observed resonance frequencies and theoretical components from the model 
described hereafter in § I4.3.1|) . 

We performed also a PCA separately for each measurement series (A and B: 9 and 5 variables, 
respectively). From series A we detected 3 important factors capturing 90.8% of the variance 
(56.9, 23.0, 11.0%, respectively). From series B we detected also 3 important factors capturing 
94.1% of the variance (54.0, 26.8, 13.3%, respectively). A 4th factor would capture only 3.6% 
more for series A and 3.4% for series B. One factor seemingly disappeared, compared with the 
PCA performed on both series. A hypothesis is that this factor is related to the hygrometric 
change between the two series. 

>From Tabled we see that all transversal and generic modes are well correlated with / actor y, 
factor2 correlates with high frequency flexural modes of both series (however notably better with 
those of series A); factory well correlates with Fl of both series (and somewhat with other low 
frequency modes: ATI, BTl and AF2), whereas factor^ correlates quite well with high frequency 
flexural modes of series B. 

3.4 Conclusions from the statistical analysis 

Fifteen modes of vibration of the clarinet reed have been observed, while previous studies investi- 
gated 4 to 5 modes only [I7J[TH]. The observed resonance frequencies are often highly correlated, 
especially those among the "transversal" modes and, to a lesser extent, those among the "flex- 
ural" modes. The nominal reed strength is surprisingly better correlated with the frequencies 
of "transversal" modes as with those of "flexural" modes. The flexural modes within the same 
series are poorly correlated. 

A principal component analysis of the resonance frequencies identifies 4 main factors, cap- 
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turing 91.2% of the variance of the sample. The data can therefore be reconstructed with 4 
uncorrelated factors only (error: RMSD = 21.8 cents, see Appendix IC.2I and [D|) . The effect of 
hygrometric change between both measurement series can seemingly be described with 1 factor 
only. 

These statistical facts offer a guidance for modeling appropriately the mechanics of the clarinet 
reed. 

4 Development and selection of a satisfactory mechanical 
model 

The determination of the mechanical parameters of a natural material is not an easy task. This 
requires to determine for each axis of the orthotropic material the value of 3 parameters (Young 
modulus, shear modulus and Poisson's ratio). These 9 parameters may exhibit viscoelastic be- 
havior, requiring theoretically for each one the fit of a viscoelastic model, such as the general 
linear solid (also called Zener model or 3-parameter model, see [2^1 HOI Ell E2] ) ■ Other multidi- 
mensional viscoelastic models could be also considered. In order to fit a wide range of frequencies 
(more than 2 decades), a 4-parameter model with fractional derivative would be required [3"2"] . 
In addition, these parameters are known to be sensitive to moisture content. Moreover, the cane 
is not homogeneous. The stiffness varies in radial direction [7] and local irregularities may be 
important, as shown by J.-M. Heinrich [33] . 

We could not consider all these aspects in the present study. Our concern is to develop a model 
with a minimal number of components, that adequately reconstructs the observed resonance 
frequencies of our reeds. We presume that these components offer an objective characterization 
of the material composing each reed. A subsequent study should examine if these components 
are correlated with some musical qualities of the reeds. This could help the reed makers to gain 
a better control on their products. 

This section is organized as follows: We first present results of FEM simulations (§ 14. ip . 
assuming an elastic and orthotropic behavior of the reed (modeled in § 14. 1 . 1[) . This helps to 
identify the modes in experiments, and allows to obtain a fit formula (§ 14. 1 .3[) . ("metamodel"), 
for computing the 11 lower resonance frequencies with respect to two parameters only, El and 
Glt, detected after a sensitivity analysis (§ I4.1.2p . This, and the previous result of 4 factors 
given by PCA, indicate that the orthotropic, elastic model is not sufficient to establish a satisfac- 
tory model with 2 degrees of freedom only. Therefore in the next subsection (S I4.2[) a viscoelastic 
model of cane is sought in order to find a fair reconstruction of the observed resonance frequencies 
for each reed and each series, after estimation of the coefficients in the model. The viscoelastic 
model f§ l4.2.2p has however too many degrees of freedom (12, for the 2 series of measurements, 
compared to the 4 factors detected by PCA), for solving adequately the inverse problem f§ !4.3J) . 
Multiple regression (see Appendix [C] for details) permits to regulate this drawback, introducing 
correlations among parameters and reducing the degrees of freedom to a number of 4. A hierar- 
chical structure can be introduced in the model, isolating the hygrometric component, bringing 
the remaining 3 components to a common basis f§ I4.3.1|) . simplifying the problem (reduced to 
only 9 "active" regression coefficients) and giving some insight in the structure of the data. The 
results are discussed in (§ I4.4p . 
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0.4087 


0.1053 


0.1835 


0.2093 


0.2235 


Glt 


0.0140 


0.2681 


0.1962 


0.1067 


0.1575 


Et 


0.0076 


0.0976 


0.0741 


0.0166 


0.0562 


Glr 


0.0438 


0.0131 


0.0257 


0.0818 


0.0341 


Er 


0.0176 


0.0120 


0.0135 


0.0662 


0.0207 


Gtr 


0.0046 


0.0092 


0.0077 


0.0215 


0.0091 


vtr 


0.0015 


0.0009 


0.0015 


0.0054 


0.0019 


vlt 


0.0018 


-0.0031 


0.0004 


0.0007 


-0.0001 


vlr 


0.0009 


0.0002 


0.0005 


0.0013 


0.0007 



Table 3: One- At-a- Time sensitivity study by FEM: Averaged ratio between relative change in 
frequency and relative change for each elastic coefficient (i.e. ±10%), sorted by decreasing order 
of magnitude, for the first 16 eigenmodes. F: flexural modes (Fl to FA), T: transversal modes 
(Tl to TA), X: generic modes (XI to X6), L: lateral modes (LI and LI; these modes were not 
observed in our study), All modes: averaged ratio over all modes. In bold: maximum absolute 
value for each coefficient of the orthotropic material: El, Et and Er: Young moduli; vlt, vlr 
and vtr'- Poisson coefficients; Glt, Glr. and Gtr' shear moduli. 



4.1 Elastic model 

4.1.1 Modeling the reed 

The clarinet reed is defined in a Cartesian axis system x, y, z (see Fig. [SJ. The origin is located in 
the bottom plane, at the tip of the reed. The material is defined as 3D orthotropic and assumed to 
be homogeneous, whose longitudinal direction L is parallel to the x axis, the tangential direction 
T parallel to the y axis and the radial direction R parallel to the z axiQ 

The dimensions in the xy plane are consistent with the measurements given by Facchinetti 
et al. [1(3] . The heel of the reed is made out of a cylinder section, diameter 34.8 mm, maximum 
thickness 3.3 mm. The shape of the reed is defined in Appendix O 

During playing, the reed has two contact surfaces with the ligature. For the present simula- 
tions, the reed is clamped in the same way than for normal playing, on two rectangular surfaces 
23 x 1 mm, spaced laterally by 5 mm, 38.2 mm from the tip of the reed, simulating the contact 
surfaces on the Vandoren Optimum ligature. However, unlike normal playing, the whole vamp 
of the reed is free to vibrate (see Figure [T]) . 

For the simulations, the "Generative Part Structural Analysis" module by Catia v. 5. 17 (Das- 
sault Technologies) is used, with mesh Octree3D, size 2 mm, absolute sag 0.1 mm, parabolic 
tetrahedrons. The generated mesh involves 5927 points, allowing both a good accuracy and a 
reasonable computing time (around 35 seconds). 

4.1.2 Sensitivity analysis of elastic coefficients 

An elastic, orthotropic and homogeneous material is defined by 9 independent coefficients: 3 
coefficients (Young modulus E, Poisson's ratio v and shear modulus G) for each axis (L: lon- 
gitudinal, T: tangential, R: radial). For selecting the most relevant parameters, we conducted 
a One- At-a- Time sensitivity analysis 34 , varying each coefficient by ±10% and computing the 
first 16 modes, based on the following reference values: El =14000 MPa, Et = Er = 480 MPa, 
vlt = v L R = v T R =0.22, Glt = 1100 MPa, G L r = G T r = 1200 MPa. The density p was set 

4 Do not confuse the morphological mode classes LI, L2, Tl, Tl, T3 and T4 with the axes L and T of the 
orthotropic material. 
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to 520 kg/m 3 , according to the estimation by Guimezanes [17]. The results are shown in Table 
131 Notice that El and Glt plays a decisive role, while Et plays a marginal role and all other 
parameters have an almost negligible influence on the resonance frequencies. As a consequence, 
the moduli El and Glt are the variables retained in the model. The approximate value of Et 
has been estimated according to the morphology of the patterns of higher order modes. This 
value is consistent with measurements given by Spatz et al. pQ. 

Notice that these results show the validity of a 2D approach, the reed being modeled as a 
thin plate. This should be used for further studies. 



4.1.3 Metamodel approximating the resonance frequencies 

The following analytic formula ("metamodel") predicts quickly and efficiently the resonance fre- 
quencies of a clamped/free clarinet reed. It was established in the following way: Frequencies 
of the first 16 modes were computed by FEM, according to a network of 92 separate pairs of 
values for E L and Glt, ranging from 8000 to 17000 MPa and 800 to 1700 MPa, respectively. 
The other elastic coefficients were held constant, according to the reference values cited above. 
For the range of simulation values, this arbitrary formula (developed by trial and error) provides 
a very good fit (generally better than ±5 cents, see Table [U]). Expected resonance frequencies / 
are first found in cents (FC) from the note F6 (1396.9 Hz), and finally in Hz: 

f(m, El, Glt) = 1396.9 x 2 FC / 1200 , where (1) 
FC i E p + a rn 2Gp + a TOi 3 E p G p + a m ^ E p + a rn ^G p , 

Ep = E L -°- 66643 and G p = Glt a762T . 

The index m is the number of the mode defined in Appendix [B] Table \§\ where the values of the 
coefficients a„ hq are given (El and Glt are expressed in MPa). 

The influence of the density is easy to predict: frequencies vary proportionally to p -1 / 2 . The 
computing cost of this metamodel is about 10 7 times lower than with FEM, largely simplifying 
the inverse problem. 



4.2 Viscoelastic model 

Equation ([T]) can be used to estimate the values of El and Glt, providing a faithful reconstruc- 
tion of the observed resonance frequencies. Theoretically these values could be computed for any 
pair of modes, after their respective observed frequencies. Unfortunately, this method gives no 
consistent results. A least squares fit is a more robust technique for such a computation. This 
leads however to systematic errors in the predicted frequencies: low-order modes are system- 
atically overestimated, while high-order modes are underestimated. This can be corrected by 
adjusting the coefficients a m fl (from Table O, but this cannot explain the bad correlation among 
flexural modes within the same series (see Table [5J ■ According to the elastic model, these corre- 
lations should be in all cases greater than 0.998. A hypothesis for resolving this contradiction is 
that the moduli are varying with the frequency in an individual way for each reed. We consider 
in this subsection a viscoelastic model, where El and Glt are frequency dependent. This leads 
to the addition of some parameters, which are to our mind more important that the other elastic 
coefficients. The fit of such a model requires many observations at different frequencies, in order 
to reduce the influence of measurements errors and of local irregularities in the structure of cane. 

Alternative hypotheses could be considered in this context, as damping, acoustic load |16j . 
local variations in stiffness or in density, local deviations in thickness, compared to the assumed 
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theoretical model. However, these hypotheses are probably unable to explain the hygrometric- 
induced individual variations we observed for each reed, thus our preference for the viscoelastic 
hypothesis. 

4.2.1 Some generalities about viscoelasticity 

It is well known, that the stiffness of natural materials like wood or cane, similarly to that of 
plastics, varies with the frequency of the applied stress and with the temperature. The material 
is stiffer at low temperature and at high frequency. At low frequency or high temperature 
the material is almost perfectly elastic and reaches the rubbery modulus. At high frequency 
or at low temperature the glassy modulus is reached; the material is almost perfectly elastic, 
also, but stiffer. At mid frequency or mid temperature, the apparent modulus (called storage 
modulus, i.e. the real part of the complex Young modulus for this frequency) is between the two 
values. For a particular frequency, called relaxation frequency, the storage modulus is exactly 
at the average of glassy and rubbery moduli. Around this frequency dissipation is maximum. 
Once the characteristic curve is known (for given temperature and different frequencies, or for 
given frequency and different temperatures), the Arrhenius equation^) offers usually an adequate 
estimate of the stiffness for any frequency and any temperature, within a quite broad range [29| . 

Considering viscoelasticity leads to complex modes with complex eigenfrequencies. Compared 
to the non-dissipative, elastic case computed by FEM, the main consequence of viscoelasticity, 
besides dissipation, is that stress and stain are not in phase. For sake of simplicity, we limit 
the computation to eigenfrequencies only, and assume that they depend on the storage moduli 
only (i.e. dissipation has a negligible influence). Having reduced the viscoelastic problem to 
an associated elastic one, the elastic solution may be used (see e.g. Ref. [IS]). In order to 
compute the resonance frequency w r after an elastic model, according to Ref. |32j . we admit 
that E ~ E'(u) r ), where E'(u) r ) is the real part of the complex modulus in the frequency domain. 
This hypothesis implies that the calculation of the eigenfrequencies from the values of the storage 
modulus is done by an iteration procedure. 



4.2.2 Choice of the viscoelastic model 

In this section, a model of cane is considered, based on the standard viscoelastic solid, also 
called Zener model or 3-parameter model (see e.g. Refs 1301 ED E2)- The scheme of the 
standard viscoelastic solid is presented on FigO with two springs E\ and E2 and a dashpot E^. 
At low frequencies, Ei and E3 have practically no effect (the rubbery modulus E\ dominates). 
At high frequencies, £3 has practically no effect (the glassy modulus E\ + E2 dominates). In 
the frequency range near E2/(2ttEq) the dissipation due to £3 is maximal and the apparent 
modulus (storage modulus) is in the mid-range. The stress a and the strain e are related by the 
constitutive equation: 

a + n& = Ei(e + r 2 e) (2) 

in which n = E3/E2 is called the relaxation time and T2 = E^{E\ + E 2 ) / (EiE 2 ) the retardation 
time. Ei is called rubbery modulus and E\ + E 2 glassy modulus. In harmonic regime, for an 

5 The shift in relaxation time is: Ln(shift) = ^ t^t — j, 1 - ^ , where E a is the activation energy, _R is the 

gas constant (8.314 J/K mol) and T and T re f the absolute temperatures in K. For instance, a shift of +10°C 
from a reference temperature of 20°C decreases the relaxation time by 16%, if E a is 13 kj/mol. 

6 This notation allows to write the parameters of the model as a vector, as required by the computations, but 
unfortunately it hides the fact that the nature of E3 (dashpot) is physically different from E\ and E2 (springs). 
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Figure 7: Schematic representation of the standard linear solid: two springs E±, E 2 and a 
dashpot £3. 



angular frequency u>, the Young modulus is complex: 



El -E 2 + jujE 3 

— = E 1 +E 2 + 

E 2 +ju)E 3 l + (u}E 3 /E 2 y 



E*(oj) = E 1+ E 2 - - : 2 — = E 1+ E 2 + , J ,\ (3) 



The second formulation separates the real part (E'{ui): storage modulus) and the imaginary part 
(E"(uj): loss modulus) of E*(u>). The storage modulus can thus be written as: 



B, <»>-* + *-:^3- ft TEW- (4) 



E'(0) = E 1 ; E'(l/ Tl ) = E!+ E 2 /2 
dE' 1 1 \ E 3 



Notice the properties: 



E'(oo)=E 1 +E 2 , 

ouj v n / 2 

This model is applied to both moduli El and Glt- For sake of simplicity the parameters are 
denoted Ei, E 2 , E3, and G\, G 2 , G3, and the storage moduli given by Equation (g]) E'(ui) and 
G'{uj), respectively. Therefore for each reed (and each series), the model requires 6 parameters 
instead of 2 (while experiments gave 4 main factors only for the whole set of results). 

>From the knowledge of the 6 parameters, the resonance frequencies can be deduced by an 
iteration procedure. For each mode the starting point of the iteration is the mean value /( ) of 
the experimental resonance frequency (see Table [T| , then the storage moduli are deduced from 
Equation then a new value /W by using Equation ([T]), etc... The convergence of the iteration 
method is fast, actually one iteration is enough. This can be understood by the fact that the 
derivative of the iterated function is small (notice that the two first rows of Table [3] correspond 
to the derivative of El and Glt with respect to frequency). If we give an arbitrary value, for 
instance = 6000, one iteration more is required for a comparable precision. In all hypotheses 
(see Appcndix lC.6[) . we used one iteration only. This procedure allows the determination of f(m) 
from the coefficients Ei, E 2 , E 3 , G\, G 2 , G 3 for a given reed and a given series. 



4.3 Inverse problem: selecting a satisfactory model 

In order to solve the inverse problem for each reed (and each series), we use a classical Mean 
Squared Deviation method, from the experimental values of the 11 resonance frequencies listed 
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Figure 8: Hypothesis H4: plot of the storage moduli E' l (uj) and G' LT (u>) [in MPa], after Equations 
Jl]to[7]), computed for the mean value of all reeds, (a): E' L series A, (b): E' L series B, (c): G' LT 
series A and (d): G' LT series B. For G" LT , the moduli are multiplied by 10. The abscissa of the 
relaxation frequencies [in Hz] are denoted by dashed lines (series A) and dotted lines (series B). 
The corresponding numerical values are listed in Tables [5] and [7] Only the portions of the curves 
between 2 and 18 kHz could be fitted adequately. The curves outside this range are purely 
hypothetical: we have no measurements. 

in Table We tested different hypotheses to establish a satisfactory model. The details can 
be found in Appendix [C] A constant value elastic model (hypothesis HI) or 2 parameter elastic 
model (hypothesis H2) are not sufficient, whereas a 12 parameter viscoelastic model (H9) con- 
ducts sometimes to non-physical results (negative rubbery modulus, for instance). For respecting 
the assumption of 4 independent parameters (as detected by PCA), correlations among the vis- 
coelastic parameters have to be introduced, as a linear combination of 4 orthogonal components. 
These 4 components are very similar to the 4 factors computed by the PCA, but, because of 
the non- linear nature of Equations ((T]and[U), a smau deviation is inevitable for optimal results. 
Factors and components are consequently strongly correlated (>0.95, see Table [2]) ■ A full linear 
model with 60 independent coefficients (H7) is not necessary, because many coefficients are very 
small. A 9 coefficient linear model (H4) is very satisfactory and we describe it hereafter. 

The RMSD (Root Mean Square Deviation: see Appendix IC.2[) is found to be 30.4 cents for 
H4, very close to 29.8 cents for H5 with 9 coefficients more. Moreover the standard deviation 
of the residuals for hypothesis H4 (and also hypothesis H5) varies very few over the different 
resonance frequencies (all around 30 cents). 

For sure, the quality of fit cannot be considered as a perfect and definitive proof that our 
models reflect the true values of the corresponding storage moduli. The influence of some missing 
parameters in the model should be examined (for instance differences in thickness between reeds, 
non constant modulus Et, non constant density p or radial variation of El). Anyways, the 
presented models reflect real mechanical differences between the reeds, very similar to those 
detected objectively by the PCA. 
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Increase in 
storage modulus 


E L 


Glt 


Series A 


Series B 


Series A 


Series B 


Mean 


31% 


16% 


38% 


36% 


Standard deviation 


5% 


6% 


10% 


9% 


Minimum 


20% 


1% 


13% 


15% 


Maximum 


44% 


32% 


52% 


49% 



Table 4: Optimal 9-coefficients model (hypothesis H4): Statistics of the frequency-dependent 
increase in storage modulus (for Series A or B and El or Glt), between the 1 st and the A th 
flexural mode: 1 - E'{2TTf F1 )/E , {2nf Fi ) or 1 - G'(2irf F i)/G , (2nf F4l ), with f F1 = 1996Hz 
and f F 4 = 16784Hz. Between both series, the correlation is 0.83 for El and 1 for Glt (as a 
consequence of the simplification of the model) 



4.3.1 Description of a satisfactory 9-coefRcients model (hypothesis H4) 

This model is designed so that no coefficient can be removed without impacting notably the 
quality of fit. It can be thought as the minimal structure allowing an adequate reconstruction 
of the observed resonance frequencies, in conjunction with the viscoelastic model Equation (T?|) 
and the metamodel Equation (JTJ) . This minimal structure makes the model more robust against 
measurements errors, but it introduces probably some bias. 

As a first step, our concern is to eliminate the influence of the moisture content and bring 
both series of measurement to a common basis. If ei[n], e2[n], e 3 [n] and e^\n] are the 4 indepen- 
dent components characterizing the mechanical properties of the reed n (these components are 
conditioned similarly to PCA as orthogonal factors: mean 0, standard deviation 1 and intercor- 
relation 0), this can be achieved by reducing the components to a number of 3 for each series s 
(ei[s,n], e2[s,n] and e3[s,n] ) in the following way: 

For series A, the components remain unmodified: 

ei [Series A, n] = ei[n] 
^[Series A, n] = e 2 [n] 

e 3 [Series A, n] = e 3 [n] (5) 
For series B, the components are reduced as follows: 

e\[SeriesB, n] = cio + ei[n] 
e 2 [SeriesB,n] = c 2 o + ^(e 2 [n] + e^n}) 

e 3 [SeriesB, n] = c 30 + e 3 [n] (6) 
The viscoelastic parameters of the model for series s and reed n can then be estimated as: 

Ei = d w + du(e 2 [s,n] + e 3 [s,n}) 
E 2 = d 20 + d 12 (2e 2 [s,n] - e 3 [s,n]) 

E3 = ^30 
G\ = g?4i(6 + ei[s, n]) 
G 2 = 2tf 41 (3 + ei[s,n] - e 3 [s,n]) 

G 3 = d eo (7) 
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This implies some other interesting relationships: 

Ex + E 2 = d w + d 20 + 3di 2 e 2 [s, n] 
Gi + G 2 = d 4 i (12 + 3ei [s,n] -2e 3 [s,n]) 
2Ei - E 2 = 2d w - d 2 o + 3di 2 e 3 [s, n] 

2G 1 -G 2 = 2d 41 (3-e 3 [s,n}) (8) 

Notice that the glassy modulus of El (i.e. E\ + E 2 ) depends linearly only on e 2 [s,n]. The 
quantities 2E\ — E 2 and 2G\ — G 2 depend linearly only on e 3 [s, n], however with opposed signs. 

The values of the 9 coefficients are: c 10 =1.011, c 20 = -2.197, c 30 = 0.82 94, di =10300, di 2 =640.5, 
d 20 =7309, d 30 =0.2822, d 41 =115.7, d 60 =0.02038. The coefficients in Equation (fTJ) are adjusted 
(in order to remove systematic errors) by adding to a mjU (from Tabled the following values, for 
m = 1 to 11: -26.27, 32.24, -50.05, 4.80, -26.87, -65.28, -48.64, 0.07, -52.98, -76.61 and -111.55 
cents. 

The change in density between both series of measurement was not precisely measured (about 
-2 to —4%). In the model, density is considered as constant. 

Fig. [5] shows an approximation of the frequency dependence of El and Glt, computed for 
the mean value of all reeds. For Series A, the storage modulus E'(uj) increases from 11700 MPa 
at 2 kHz (Fl) to 17200 MPa at 16.8 kHz (F4), while for Series B, this increase goes from 11100 
to 13300 MPa. Under "normal" hygrometry the reed bends therefore in a notably viscoelastic 
manner, whereas the ultra-dry reed bends in an more elastic manner. For Glt, the reed is 
generally always notably viscoelastic, according to our model. The corresponding values are: 
783 to 1290 MPa for Series A, and 896 to 1436 MPa for Series B. Corresponding statistics are 
displayed on Table 01 Drying seems to increase Glt and decrease El (except around 2300 Hz), 
explaining the good correlation between the variables AF1 and BF1. 

This simplified model permit interesting conclusions about the structure of our data: 

• The component e\ is related exclusively to Glt] 

• e 2 is related exclusively to El; 

• e 3 increase proportionally the rubbery modulus of El (and not its glassy modulus), and 
decrease proportionally the glassy modulus of Glt (and not its rubbery modulus); it 
regulates therefore the viscous component common to El and GltQ; 

• e 4 takes into account the variation in moisture content between series A and B. 

Notice that the product of the rubbery moduli for the series A, E\G\ is correlated at 72% 
with the nominal reed strength. Individually, these moduli correlate only at 50% and 55% with 
the nominal reed strength, respectively. 

4.4 Discussion 

Our viscoelastic model is able to explain partially the bad correlations observed between flexural 
modes. In Table[5J we compare the correlations among observed resonance frequencies and among 
computed modal frequencies, after viscoelastic model, hypothesis H4 and H5. It seems that 
additional hypotheses (such as an irregular thickness) should perhaps be considered for improving 

7 The parameter E2 (respectively G2) determine the influence of the Maxwell arm in the Zener model, since 
-E3 (respectively G3) is constant in our simplified model; if E2 = the model is perfectly elastic and the viscous 
component disappears 
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Model 


Mode 


AF1 


AF2 


AF3 


BFl 


BF2 


BF3 




F2 


0.73 






0.63 






observations 


F3 


0.34 


0.75 




0.43 


0.84 






FA 


0.28 


0.65 


0.74 


0.00 


0.43 


0.45 




F2 


0.93 






0.79 






H4 


F3 


0.69 


0.90 




0.62 


0.97 






FA 


0.56 


0.81 


0.98 


0.36 


0.82 


0.91 




F2 


0.90 






0.69 






H5 


F3 


0.57 


0.87 




0.52 


0.98 






FA 


0.43 


0.77 


0.98 


0.26 


0.84 


0.92 



Table 5: Correlations of resonance frequencies between flexural modes within the same series, 
after observations and viscoelastic models H4 and H5. AF1, AF2, AF3, BFl, BF2 and BF3: 
correlations with modal frequencies Fl, F2 and F3 within series A or series B. Lines: corre- 
sponding flexural mode between which the correlations are computed. 

the model. However, we should remember that the determination of resonance frequencies are 
attached with some uncertainties, especially for the modes F3 and T3. 

In order to clarify this question, let us examine if the residuals (observed resonance frequencies 
minus computed modal frequencies with hypothesis Hfjf|) contains some pertinent information. 
A PCA shows that perhaps 2 residual factors contains some interesting information (explaining 
30% and 18% of the residual variance). The first residual factor is correlated with AT3 (0.86), 
BT2 (0.74), AT2 (0.73) and AX2 (0.66). These modes depend all strongly from G LT . An 
adjustment of the coefficient a m , 2 from the metamodel for the transverse modes could probably 
cancel this systematic bias (remember that the coefficients a m>q are computed from a theoretical 
model, which is probably also biased). Indeed, an increase of 14, 21, 16 and 11% of this coefficient 
affecting the modes Tl, T2, T3 and TA drop the RMSD from 29.8 to 28.5 cents. The second 
residual factor is correlated with AF1 (0.36), AF2 (0.29), AT2 (-0.29), BFl (0.28) and BT1 
(—0.23). This reveals probably a competition between flexural and transversal modes when fitting 
the model. The bias comes probably from the coefficient a mj i, regulating the linear dependance 
to El in the metamodel. Adjusting the coefficients a m) i for all modes and the coefficients o mj 2 
for all transverse and generic modes drop the RMSD down to 26.2 cents, reaching practically 
the size of the measurement steps (25 cents). The adjustments of a m .i are very small for the 
flexural modes Fl to FA: -6, 1, 6 and 1%. 

This shows that the most important bias depends linearly from the 2 most important pa- 
rameters (El and Glt) of our FEM computations. No supplementary parameter is required 
until this bias is removed (theoretically down to a RMSD of 21.8 cents, according to H10). For 
hypothesis H4, the same linear adjustment of the metamodel let drop the RMSD from 30.4 
down to 27.0 cents. 

A comparison with results by other authors is difficult or even quite impossible, because of 
the disparate structure of the measurements. Such a comparison requires the reconstruction of 
the measurement data (when possible), the fit of a viscoelastic model or the extrapolation of the 
values, in order to reach the frequency and temperature range of our measurements. The validity 
of such a highly speculative task is questionable. For the storage modulus E' L , the reconstructed 
values from other authors fall all in the range around the mean of our measurements ±3 times 
the standard deviation, however most of the time in the lower range. This shows probably that 
the selected value for the density p was somewhat too high. No representative statistics are 

8 Because H5 is probably less biased than H4 with its minimal structure 
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Model 


Storage Young modulus E' l {uj) at about 20°C 


E x [MPa] 


E 2 [MPa] 


fr [Hz] 


E' L at 4 kHz 


Hypothesis H4, series A 


10300, 
SD 906 


7309, 
SD 1432 


4123, 
SD 808 


13781, 
SD 897 


Hypothesis H5, series A 


9377, 
SD 844 


8336, 
SD 1552 


4027, 
SD 750 


13449, 
SD 839 


Hypothesis H4, series B 


9423, 
SD 784 


3964, 
SD 1109 


2236, 
SD 626 


12338, 
SD 885 


Hypothesis H5, series B 


7844, 
SD 753 


5459, 
SD 1217 


1947, 
SD 434 


12168, 
SD 882 



Table 6: Summary of our results for hypotheses H4 and H5 about the viscoelastic behavior of 
the longitudinal Young modulus El in cane. E\, E 2 and f r (relaxation frequency): parameters 
from Zener model. E' L at 4 kHz : storage modulus at 4 kHz [in MPa]. SD: standard deviation 
(the value preceding SD is the average). The model is valid only between 2 and 18 kHz. 



Model 


Storage Shear modulus G' lt (lu) at about 20°C 


Gi [MPa] 


G 2 [MPa] 


U [Hz] 


G' LT at 4 kHz 


Hypothesis H4, series A 


694, SD 116 


694, SD 327 


5420, 
SD 2555 


926, SD 102 


Hypothesis H5, series A 


752, SD 119 


628, SD 328 


6310, 
SD 3296 


924, SD 105 


Hypothesis H4, series B 


811, SD 116 


736, SD 327 


5749, 
SD 2555 


1042, SD 99 


Hypothesis H5, series B 


774, SD 119 


769, SD 328 


5622, 
SD 2401 


1022, SD 103 



Table 7: Summary of our results for hypotheses H4 and H5 about the viscoelastic behavior of the 
shear modulus in longitudinal / tangential plane Glt in cane. G' LT at 4 kHz: storage modulus 
at 4 kHz. Same structure as Table O The model is valid only between 2 and 18 kHz. 



available by the other authors. 

The most important disagreement, compared with our model (hypotheses H4 and H5), is 
about the relaxation frequency. The explanation is probably because the studied frequency 
range was not same. For Glt, we found no viscoelastic measurements by other authors. Our 
results are summarized in Tables and [7] For E' L , between hypotheses H4 and H5, f r and E' L at 
4 kHz agree well, whereas E\ and E 2 diverge by about 1-2 SD. This divergence comes because 
the observed frequency range was not broad enough. For the shear modulus, G\ and G 2 are in 
good agreement for both hypotheses (and consequently f r and G' L at 4 kHz also). 

Our model is valid only for "ambient dry" reeds (since the ultra-dry conditioning was not 
controlled), in a frequency range which should not exceed one decade. We checked that a 
fractional derivative model after Gaul et al. |32j is not necessary in our narrow frequency range. 
Such models are however really efficient to cover a broad frequency range. For instance, the 
data by Lord for dry material QJ] could be fitted very well (E = 8108 MPa, E x = 2964 MPa, 
p = 0.298 MPa-s Q , a = 0.546). Notice that the order of the derivative (a) is 1 in our viscoelastic 
model. 
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5 Conclusion 



The main results from the statistical analysis are discussed in WS.4\ >From these results it was 
possible to establish a satisfactory numerical model, selecting the most important parameters 
describing the mechanical behavior of a reed. 

A sensitivity analysis by FEM calculation assuming an orthotropic, elastic material has been 
conducted, and showed that the longitudinal Young modulus El and the longitudinal / trans- 
verse shear modulus Glt play a leading role. An efficient metamodel could be developed to 
approximate the resonance frequencies of a clarinet reed from these two moduli. In order to 
improve the model, we consider the hypothesis of frequency dependent moduli, introducing ad- 
ditional degrees of freedom. A reconstruction of the observed resonance frequencies can then 
be achieved with a good accuracy, estimating for each reed only 4 components, from which the 
parameters of a viscoelastic model are computed as a linear combination. The selected model 
(after hypothesis H4) is probably slightly biased, but it is more robust against measurement 
errors than more refined models. 

Table [H shows that these components are highly correlated to the factors computed by PC A 
(0.96 to 0.98). A multiple regression confirmed the data structure detected with the viscoelastic 
model. 

The proposed method allows the determination of 3 mechanical parameters characterizing 
the material composing each reed, with a single series of measure, using Equations JTJ U] and 
[7]). The reed should be conditioned with a relative humidity corresponding to the one ensured 
by the hermetically sealed package by Vandoren (about 55%). The forth parameter cannot be 
determined in a reproducible way, since the exposition of the reeds to the ultra-dry air of the 
optical laboratory was not controlled. The same protocol and the same viscoelastic model can 
be used for other kinds of single reeds (bass clarinet, saxophone). Only the coefficients of Table 
[5] have to be recomputed after a FEM simulation of the corresponding reed shape. 

Despite the fact that the eigenmodes of higher order probably play no important role in the 
acoustics of the clarinet, the present study shows that they reveal the inner structure of the 
material building the tip of the reed, so a new step could be done for an objective mechanical 
characterization of the clarinet reed. 
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APPENDICES 



A Denning the shape of the reed 

The thickness of the vamp at point (x, y) is interpolated as follows: first, we interpolate 3 points 
at y=0, 4 and 6mm, with 3 cubic splines, after Table [5] These three points (sq(x), s\(x) and 
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x [mm] 


So [mm] 
y=0 mm 


si [mm] 
y=4mm 


So fmml 
y=6mm 





0.074 


0.080 


0.042 


5 


0.343 


0.293 


0.197 


10 


0.648 


0.542 


0.377 


15 


1.047 


0.847 


0.571 


20 


1.451 


1.135 


0.745 


25 


1.926 


1.527 


1.078 


30 


2.540 


2.084 


1.589 


35 


3.351 


2.817 


2.256 



Table 8: Network of points for interpolating the thickness of the vamp. See explanations in the 
text. 



S2(x)) define a biquadratic polynomial: 

vamp(x,y) = p (x) + pi(x) y 2 +p 2 (x)y 4 

withpo(a;) = So{x), Pi (x) = (— 65 s (a;) + 81 s\(x) — 16 S2(x))/720 &ndp2(x) — (5s (x) — 9sl(x) + 
4s2(2;))/2880, allowing an interpolation on the y axis. The network of points above was estimated 
using a least squares fit, based on a network of 12 x 24 thickness measurements, achieved with a 
dial indicator and a coordinates-measuring table (estimated accuracy: ±5^to in z, ±50/ito in x 
and y). We measured twenty reeds and select a particularly symmetrical one as reference. This 
method allows the reconstruction of the measurement network with an accuracy of ±10/im. 
The thickness of the heel is defined to be : 



heel(y) = -14.1 + ^17 A 2 - y 2 
The contour of the reed in the xy-plane is defined by: 



contour(x) = < 



fo 

v/(24.4 -x)x 
4.08044 - 

263 11 



V-5.31 + 6.8x - a; 2 



V 40 900 

The thickness of the reed at point (x, y) is defined by: 



x < or x > 67.5 
x < 1.13196 
x < 2.94661 
x < 67.5 



thickness(x, y) 



min[heel(y),vamp(x,y)] Abs(y) < contour (x) 
otherwise 



B Coefficients of the metamodel 
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mode 


m 


O"m,0 






am, 3 






5- 


s+ 


Fl 


1 


2334.56 


-1165877 


0.2763 


-21.99 


145883795 


-0.000324 


-2 


4 


Tl 


2 


1481.70 


-642027 


5.0725 


569.69 


83451520 


-0.005886 


-3 


3 


F2 


3 


3651.79 


-1060625 


0.7082 


-123.77 


130061359 


-0.000786 


-2 


5 


T2 


4 


2403.82 


-493591 


3.8130 


501.70 


64716271 


-0.003741 


-9 


5 


XI 


5 


3153.73 


-822597 


3.1425 


604.29 


105032670 


-0.003785 


-3 


4 


F3 


6 


4669.87 


-1009348 


0.7623 


-101.18 


122741172 


-0.000877 


-3 


6 


T3 


7 


3015.77 


-275531 


3.5285 


250.00 


29933670 


-0.003016 


-2 


1 


XI 


8 


3874.44 


-633268 


2.0921 


589.25 


83462866 


-0.001958 


-9 


6 


X3 


9 


4381.85 


-926543 


1.9907 


730.31 


122925338 


-0.002593 


-3 


3 


T4 


10 


2659.08 


588457 


6.6088 


-1269.38 


-120493740 


-0.005269 


-19 


26 


FA 


11 


6450.01 


-1689363 


-2.5544 


1282.75 


247738162 


0.001631 


-32 


42 



Table 9: Coefficients of the metamodel, Equation ([T]). The maximum negative and positive 
deviations of the model (compared to the values calculated by FEM) are given by S— and 5+ 
[in cents]. The mode LI (6th mode) has been deleted from the Table, since we didn't observe it. 
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Indice 


from 


to 


numbering 


n 


1 


N=55 


reeds 


s 


1 


S=2 


series of measurements A (s=l) and B (s=2) 


m 


1 


M=ll 


modes 


q 





Q=5 


coefficients for Equation ([T]) 


i 


1 


1=2 


moduli E L (i=l) or G LT (i=2) 


j 


1 


J=3 


viscoelastic coefficients 


k 





K=2,3 or 4 


components (factors) 



Table 10: List of indices 



C Development and selection of a simplified viscoelastic 
model 

In this section we describe how the proposed model was developed and selected. Alternate 
options are presented. 

C.l Data structure 

We need a specific notation for denoting our complicated multivariate data structure as arrays, 
after a list of indices (see Table [TU|) . The 11 modes are defined after Table [5] (from 1 to 11: Fl, 
Tl, F2, T2, XI, T3, F3, X2, X3, T4, FA). We define a variable v ni8>i)j holding all coefficients 
of our viscoelastic model: 

v n ,s,i,j = Ej and u n ,s,2,j = Gj for reed n and series s. (9) 

The array r njSjlrl holds the reconstructed resonance frequencies computed for all reed, series 
and modes, with the coefficients v = {u n ,s,i,j} our viscoelastic model (see jj4.2.2|) . 

C.2 Mean Squared Deviation 

As function to minimize, we define the Mean Squared Deviation MSD (also called Mean Squared 
Error) between reconstructed and measured resonance frequencies On,s,m" 

^ N S M 

MSD = nsm EEE - r ™-») 2 ( 10 ) 

n— 1 s— 1 rn—1 

With Equation p0|) . the components of array v can be fitted by any appropriate algorithm for 
minimizing a multivariate function. All available measurements can be utilized for fitting the 
models^. Missing observations o mtSjn are then eliminated while computing MSD. 

Our model allows a quite good reconstruction of the resonance frequencies, with a V MSD = 
RM SD (Root Mean Square Deviation) smaller than 20 cents (it is very small for lower modes, 
and always smaller than 25 cents). This correspond to the hypothesis named H9 (all hypotheses 
HI to Hll are presented and commented in § IC.6[) . but the values of the coefficients are not 
always plausible physically. This problem comes out because the observed resonance frequencies 
are far from 0, so the rubbery modulus E\ cannot be estimated precisely. The algorithm tends 
to produces sometimes unlikely values. We can correct this problem by setting that the damping 
coefficients E3 and G3 are constant, independent of the reeds and eventually independent of the 

9 This was not the case with PCA. 
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series (typical values: E3 = 0.28 and G3 = 0.02). It remains 4 coefficients for each series and 
reed, E\, E 2 and G\, G2 (i.e. 8 coefficients for each reed). This ensures generally that these 
coefficients fall in a plausible range, when fitting the model. 

However, our model has still too much degrees of freedom, compared to the 4 or 5 factors 
detected by PCA. In such a situation, measurement errors and local irregularities in the structure 
of cane could have a marked influence by fitting the model. In § IC.3[ we examine how to regulate 
this drawback by multiple regression. 



C.3 Estimating the parameters of the viscoelastic model by multiple 
regression 

Multivariate linear regression consists in projecting linearly a n dimensional space on a onedi- 
mensional space. The generic equation is 



N 



y = dp + } a n x n = 2J a n x n , if x a = 1 



n=0 



Equation (JTTJ) can be generalized for multiple regression as: 

y = Ax 



(11) 



(12) 



where x and y are column vectors and A a matrix. In what follows, in order to use conventional 
vectors and matrices (with respectively one and two dimensions), we use the following notation: 
Vj[n, s, i] — v ntS _i t j is the jth component of the column vector v[n, s, i]. 

In our case, we have no prior knowledge about the relationships among the parameters of our 
model. We assume that each parameter in the model can be computed as a linear combination of 
some unknown independent components by multiple regression. The multiple regression formula 
can be written as: 



Vj [n, s, i\ 



K 

k=a 



jk[s,i\ ■ e-k[n\ 



(13) 



In conventional vectors and matrices notation, Equation (|13|) reads: 

v[n, s, i] = M[s, i] ■ e[n] 

where e[n] is the vector of the orthogonal components for each observed reed n, M[s,i] is 
the regression matrix, independent of the reed number, depending on the series and the kind of 
modulus (El or Glt)- v[n, s, i] is the vector of the coefficients of the viscoelastic model. 

We have only to choose an arbitrary number of components, for instance K = 4, referring to 
our PCA. We introduce arbitrary constraints in order to obtain more comparable results: over 
the different reeds, we state that the components must be orthogonally normalized (mean 0, 
standard deviation 1 and intercorrelation 0) . The matrix of components e has consequently to 
satisfy: 



/ N 



e ■ e = 







N-l 






V 






N-l 









N-l 




\ 





N-l J 



(14) 
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Each individual column vector e[n] from this matrix is written as follows (for the component 
eo[n] = 1 : see Equation (TTT]) '): 

e[n] = ( 1 ei[n] e 2 [n] e 3 [n] e 4 [n] ) T . 

The components of M are fitted by minimizing Equation (jlUp. As starting value for the 
orthogonal components we set : Ck [n] = f actor k [n] , where f actor k [n] are the factors computed by 
PC A (see ^3.3[) . After a first estimation of M, it is possible to release the approximation about e: 
all components and all coefficients in the matrices can be fitted by the fitting procedure. However 
the number of variables to fit is probably much higher than allowed by the most algorithms of 
function minimizing. The fitting procedure has to be carried "by hand" with subsets of variables. 
The procedure we used is described below. 

For K — 4 (the dimension of vector e[n] being 5), and some supplementary choices, the model 
works very well. 

C.4 Empirical simplification of a 4-parameter model 

We observed that the effect of hygrometric changes between both measurements series can be 
taken into account with one parameter only, practically without drop in quality of fit. This 
effect can be isolated on component e^n] and the remaining components can be transformed 
linearly, so the further computations can be achieved from a common basis. Equation (|13[) is 
then simplified in: 

M[s,i] =M[i\-M[s] , (15) 

The reasoning considers the situation where the two series of observations are independent, with 
a number of components in vector e[s,n] = M[s] ■ e[n) reduced to K = K — 1 = 3. As a 
consequence the number of rows of matrix M[s] is 4, as well as the number of columns of M[i]. 
This approach, allowing to separate hygrometry effects, offers a comfortable way to test different 
hypotheses, without changing the structure of the computation, by setting some coefficients in the 
matrices at some arbitrary values or by introducing some linear relationship between coefficients. 
We get fewer "active" coefficients to fit in the model: the fitting procedure is much faster and 
this raises the probability to find the best possible fit. 

We tried to minimize the number of coefficients different from zero in the matrices, without 
substantial drop in quality of fit. The model was fitted "by hand", using the solver of Excel 
(Microsoft Office). We found empirically that a quite sparse setup still gives a good filF°l: 

(rfio d 12 di 3 \ I d 40 du \ 

d 20 d 22 d 23 ; M[GLT] = d 50 d 5 i d 53 (16) 
d 30 / V d 60 / 

10 The fitting process was realized by performing four kinds of procedures, in an arbitrary order : 

1. Fit the "active" coefficients in model, without adjusting the coefficients a m ,o 

2. Fit the "active" coefficients in model and adjust the coefficients a m ,o (or fit all coefficients: a mj9 ). 

3. Fit the e n ^ components (individually for each reed, irrespective of Equation H14I 0, then then normalize and 
orthogonalize components (in order to respect Equation 1141 ). finally rotate components (for improving 
the fit). 

4. Eliminate some "active" coefficients (set them as 0, 1 or some other constant value; set some arbitrary 
linear dependence from other coefficients) 



2.x 



M[SeriesA] = 



M[SeriesB] = 



/ 1 











^ 


ClO 


1 











C20 





C22 





C24 


V c 30 








1 


) 



(17) 



This corresponds to the hypothesis named H5. Furthermore some coefficients may be propor- 
tional to others, without noticeable drop in quality of fit. This diminishes the number of "active" 
coefficients in the different matrices from 18 to 9 (hypothesis H4): 



c?22 = 2e?i2 ; di3 = <ii2 = — 1^23; 

cLto = d 50 = 6d 41 ; d 5 i = 2d 41 = -d 53 ; (18) 
C22 = c 2 4 = 1/2. 

C.5 Adjusting coefficients for removing systematic errors 

After fitting the different models, we observed some systematic deviations in the resonance 
frequencies between model and observations. This error has probably two different origins: 
an inevitable inaccuracy in the FEM computation (and in our metamodel) and an error for 
parameters not included in the model. A straightforward way to minimize the residuals is to 
adjust the coefficients a„ lj0 in Equation Fitting all coefficients in Equation JT]) is doubtless 

a more questionable way to reduce this error (hypotheses H8 and H9). This can reduce the mean 
deviation between model and observations, but greatly increases the number of coefficients in 
the model (see however the discussion jj4.4p . 



C.6 Testing different hypotheses 

We tested different hypotheses with our model, in order to select a particularly efficient model. 
Some of them are summarized in the Table 111! 

Hypotheses HI to H3 (elastic model) present a poor fit; the adjustments for coefficients a m ,o 
are large, compensating partially the missing viscoelastic components. All viscoelastic models 
are notably better and exhibit smaller adjustments for a m .n- Frequency dependence for El and 
Glt seems evident. Hypothesis H4 shows a good accuracy, with only 9 fitted coefficients (and 11 
adjustments). Augmenting the number of coefficients up to 60 brings only a marginal contribution 
(H5 to H7). Adjusting the other coefficients of Equation ([TJ) in H8 and H9 improve the model, 
especially for the high modes. Our FEM computations (and consequently our metamodel) are 
probably attached with systematic errors in this frequency range. The influence of Et should 
eventually be considered. Between H8 and H9, the total number of components (NxK) increases 
from 220 to 660. The adjustments within morphological classes are related: "flexural" modes 
shows systematically lower values than neighboring "transversal" modes. 



C.6.1 Backward random validation 

Following a suggestion by a referee, we applied a principal component analysis to simulated data 
computed after hypothesis H5: we assigned randomly a value for the 4 components and 55 reeds, 
following a normal distribution. We repeated this operation ten times. As expected, the PCA 
detected 4 factors capturing 91.2% (Standard Deviation 1.3%) of the variance of the simulated 

11 These coefficients can be fitted through the fitting procedure (some constraints have however to be introduced, 
to avoid an important deviation from their theoretical values) or merely adjusted, so that the total averaged 
deviation for each mode and both series is 0. 
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Hypothesis 


K 


model 


#ActiveCoef 


#OtherCoef 


RMSD 


HI 





elastic 


2 


«m,0 — > 11 


76.2 


H2 


2 


elastic 


4 


a m ,o -4 11 


54.8 


H3 


3 


elastic 


7 


a m ,o -4 11 


43.9 


H4 


4 


viscoelastic 


9 


Om,0 — > 11 


30.4 


H5 


4 


viscoelastic 


18 


Om,0 — > 11 


29.8 


H6 


4 


viscoelastic 


44 


flm,0 — > 11 


29.1 


H7 


4 


viscoelastic 


60 


flm,0 - > 11 


28.6 


H8 


4 


viscoelastic 


60 


a m , g -4 66 


23.2 


H9 


12 


viscoelastic 




a m , q -4 66 


19.8 


H10 


4 


regression 




W s , m , fe -4 110 


21.8 


Hll 


4 


regression 


4 


W m ; k 44 


24.7 



Table 11: Synthesis of some hypotheses tested with our model. K: Number of components 
for each reed n; model: elastic, viscoelastic or multiple regression (with the elastic model, £x 
and Glt are independent from frequency; the elastic model is computed after the viscoelastic 
one by setting very small values for the coefficients affecting E2, E3, G2 and G3, in order to 
avoid division by 0; multiple regression: see Equations (fTTH and |2"0|) ) : #ActiveCoef: number 
of "active" coefficients in the matrices M or M and M (fitted through the fitting procedure); 
#OtherCoef: number of other coefficients estimated in model (a mj o: adjusted so that the 
mean error for each mode m and both series is zero, W Stmi k- computed analytically otherwise: 
fitted through the fitting procedure); RMSD [in cents] after Equation (|1U|) : this is a measure of 
goodness of fit (remember that the vibration patterns of the reeds were observed in steps by 25 
cents). The results are commented in the text. 

data (mean: 42.6, 20.9, 16.2 and 11.5% for each factor). This seems compatible with the analysis 
performed on the observed frequencies f§ !3.3p : 4 factors: 91.2% of the variance (53.6, 21.4, 10.8 
and 5.4% for each factor). 

D Reconstructing observed resonance frequencies by mul- 
tiple regression 

The scheme of our viscoelastic model is: 

(vector e) -4 {multiple regression Equation (jl3p| -4 (viscoelastic coefficients v) — > {viscoelastic 
model Equation ([4])} — > (moduli E' L and G' LT ) — > {metamodel Equation ([1])} -4 (reconstructed 
resonance frequencies r). 

It has a very interesting property: the same model of cane can be used for any kind of reeds 
(for instance bass clarinet or saxophone) or for any other boundary conditions. Only Equation 
dTJ) has to be changed (or at least, the coefficients a mA have to be recomputed). 

Within our particular setup, the viscoelastic model is however not required for reconstructing 
the observed resonance frequencies: PCA is theoretically the optimal linear scheme, in terms of 
least mean square error, for compressing a set of high dimensional vectors into a set of lower di- 
mensional vectors and then reconstructing the original set by multiple regression. The shortened 
scheme is merely: 

(vector factor) — > {multiple regression} -4 (reconstructed resonance frequencies r) . 

Let us examine this option. For this purpose we use the array factor n ,k computed in § 13. 3[ 
holding our 4 principal components (factors). As before f $C.3[) . we set factor n ,o = 1. The array 
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of reconstructed resonance frequencies r njSjlri can be computed by multiple regression using an 
array of matrices W s ^ m ^\ 

r[n,s] = W[s] ■ factor[n] (19) 

As before f 9C.4p . we have also the option to reduce the dimensionality from K to K, using the 
previously defined array of matrices M[s] and then use a unique matrix W m ^ to operate the 
multiple regression: 

r[n, s] = W ■ M[s] ■ factor[n] (20) 

We call these 2 options: hypotheses H10 and Hll. For Hll, we performed a small orthogonal 
rotation of the factors to concentrate the information about the hygrometric material properties 
in f actor for a better fit. The results are summarized in Table [TT1 

Multiple regression is an accurate way to retrieve our measurements, comparable to H8 and 
H9. With only 48 coefficients, Hll is very efficient, even better for "transversal modes" as H10. 
The results of the regressive model are more difficult to interpret as those of the viscoelastic 
model. As e nj 4 before, factor n ^ serves uniquely to adjust factor n ^ relatively to series B. As 
expected, factor n ^\ influences mainly the "transversal modes" and factor^ the "flexural modes". 
Their respective coefficients in the matrix W reflect this antinomy. "Transversal" and "flexural 
modes" are concerned by factor n ^ in a quite similar way, but the slope is not same. 
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